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Thomas-Fermi Statistical Models of Finite Quark Matter[2] 
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I introduce and discuss models of finite quark matter using the formalism of the Thomas- 
Fermi statistical model. Similar to bag models, a vacuum energy density term is introduced 
to model long distance confinement, but the model produces bound states from the residual 

^\ color Coulomb attraction even in the absence of such a term. I discuss three baryonic 

O 

f~^ applications: an equal mass nonrelativistic model with and without volume pressure, the 

ultra-relativistic limit confined by volume pressure, and a color-fiavor locking massless model. 

C^ These models may be extended to multi-meson and other mixed hadronic states. Hopefully, 

| .^ it can help lead to a better understanding of the phenomenology of high multi-quark states 

Cn in preparation for more detailed lattice QCD calculations. 
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I. INTRODUCTION 

A long standing question in the theory of quark matter asks: where are the mesons, baryons or 
mixed states with more than the usual two or three quarks? It is likely from the theory, quantum 
chromodynamics (QCD), that high multi-quark states exist. A number of "exotic" hadronic objects 
have been postulated. There has been speculation concerning hypothetical strange quark statesfU 
O [H] for many years, partly motivated by astrophysical considerations involving the so-called GKZ 
cutoff for suppression of high-energy cosmic rays|i4 . It hasn't been clear until recently that the 
GZK cutoff mechanism was producing the expected suppression of events at energies above ~ 10'^'^ 
GeV^Sj. 

Experimental evidence for a four quark one antiquark state has been much discussed[6] since 
the initial report in 2003[7^, and there is a very recent report of a double meson four quark state 
from the Belle collaboration [8] . The best way to investigate such theoretical questions is with 
the set of methods known as lattice QCD. Although it is in the initial stages of being applied to 
larger nuclear systems lattice QCD is, however, still limited to small volumes and quark numbers. 
There is a need for models which can help lead expensive lattice QCD calculations in the right 
direction in the search for high quark states. One such treatment is the MIT bag model, which has 
been applied to the problem of modeling hadronic bound states with many quarks[2l HI |9]. The 
Nambu-Jona-Lasinio model has also been used in such analyses |10j. 

I would like to point out here that the Thomas-Fermi (TF) treatment of particles as a gas at zero 
temperature [m [12] provides a versatile and adaptable set of models which can be used to study 
such questions. It has been quite successful in treating atomic systems, viewing those as a bound 
collection of interacting gas particles, even for low atomic numbers. One of its attractive features is 
the exact treatment it affords to the Coulomb interactions. It can be understood as a semi-classical 
approach which correctly incorporates fermion statistics for a large number of particles but is not 
fully quantum mechanical. It has traditionally been used to model ground state atomic systems 
(it has also been applied to nuclear and condensed matter problems), but can also be applied to 
quark systems. We will find that this is a natural extension of the method. 

The following model combines the TF teatment of quarks with a bag model assumption to model 
confinement. It removes the central nuclear force present in the traditional TF atomic model and 
replaces it with a collective residual color Coulombic interaction which provides binding even in the 
absence of a volume term. Although it is related to bag models, the TF quark model has a different 
starting point. Bag models represent exact solutions to free fermions placed in a spherical cavity. 



In the TF model, one deals not with individual modes but with a single, collective, spherically 
symmetric density of states which represents a gas of 3^ ("A" is baryon number) interacting 
fermions at zero temperature. The spherical symmetry assumption is natural for a gas of particles in 
their ground states and is a prosteriori justified by bag model fits, with reasonable phenomenology. 
The quarks interact via their residual color Coulombic interactions. (The word "residual" is meant 
to denote the non-primary, non-confining color Coulombic part of the interaction. This interaction 
scales like ~ 1/A in a color singlet, and is thus residual in a large A sense also.) In contrast, bag 
models assume no interactions to lowest order and the color interactions are put in perturbatively. 
In other words, the advantage for TF models is the nonperturbative inclusion of the Coulombic 
interactions; the disadvantage that one is no longer working with exact solutions of the Dirac 
equation but assumes a statistical treatment. Nevertheless, one would expect that the TF quark 
model would become increasingly accurate as the number of constituents is increased, that is, as a 
statistical treatment becomes more justified. Note we are talking only about ground states in this 
context, as the TF model is incapable of describing excited atomic states. The main usefulness of 
the model will be in seeing systematic trends and effects in ground states as the parameters of the 
model are varied and helping to identify likely stable bound states. 

The different starting point for the TF quark models results in a very different consequence 
relative to bag models: the residual color Coulombic interactions produce a bound state in the 
massive version of the model. The attractive collective potential is shown in Fig. [3] below for the 
nonrelativistic single flavor case. However, these Coulombic bound states are very different from 
normal hadronic states. Their energies scale like A^'^ rather than A. The natural interpretation 
of these systems is that they are bound states in a zero temperature version of the quark-gluon 
("deconfined") phase of QCD. 

In order to describe normal hadrons, it is also natural to assume a phase transition to a confined 
state. This is implemented by the introduction of the bag energy, BV {V is the volume), which 
produces an external pressure and results in a surface discontinuity in the particle density function. 
Then, as a result of the total energy minimization, one finds these energies are proportional to A 



for A » 1. The nonrelativistic model has a total energy as a function of A given by Eq.(71 ). This 
equation is similar to the liquid-drop nuclear model binding energy equation. Deviations from the 
~ A behavior should tell us where to expect greater or lesser binding. The (massless) relativistic 



model's energy equation, Eq.(95), also scales like A for A » 1. 

In the following section I will first work out the Coulombic couplings associated with a color 
singlet gas of quarks whose charges are not fixed. The interactions in the TF quark model are 



modeled on the basis of classical QCD, as explained for example in Ref.|13]. As mentioned above, 
the effective couplings, which incorporate color interaction probability factors deduced from the 
classical theory, decrease like 1/^4 for large values of A. The energy functional will then be for- 
mulated in Section III where I will show that the energy equations can be formally integrated and 
related to endpoint values and derivatives of the TF function. I consider a nonrelativistic form 
of the model with Nj mass degenerate quarks in Section IV. A stable state is formed, which, as 
discussed above, is interpreted as a hadronic bound state in the zero temperature quark-gluon 
phase. In Section V I examine the effects of turning on the external pressure, and we will see 
that it produces a surface discontinuity in the TF functions. An explicit expression for the system 
energy will be derived. I will briefly consider the general form of the TF relativistic quark equation 
in Section VI before solving the specific case of Nf massless quarks. This system is subjected to 
external pressure and an expression for the system energy is derived. To show the versatility of 
the model and point to future applications, I also consider aspects of a color-flavor locking type 
interaction in Section VII, which includes a Cooper pairing term and massive gluons. I close with 
a summary and some general comments about directions of future work in Section VIII. 

During final editing of this paper, the author became aware of the paper |14J, which uses an 
approximate TF expression for the fermion partition function to evaluate the equation of state for 
a gas of massless quarks interacting with SU{2) or SU{2>) color gluons. I will comment more on 
this paper and its relationship with the present model in Section VI. See also Refs.fTS] I16|. which 
applies a TF approximation to the evaluation of the lattice QCD fermion partition function. 

II. RESIDUAL COULOMBIC COUPLINGS 

Quark-quark and quark-antiquark SU{2>) color interactions at the classical level can be un- 
derstood from their "charges" in isospin/hyper charge space, usually designed as the "3" and "8" 
axes|13]. Quarks can be understood to have strong interaction charges with magnitude g5^, where 
g is the strong interaction coupling constant, located at the corners of an equilateral triangle with 
vertices at (—1, -yn)g, (1, ~7n)9^ ^"^^ (0) ~~7n)9 iii (3,8) space notation. Likewise, antiquarks can be 
understood to exist at the points (—1, — 7?)5) (1) — 71 )5' ^^"^ (^' ^^9' ^"^ *^^^ section I will treat 
only multi-baryon states and will assume that all bound states have zero net color. Introducing 
the total charge as 

^ = ^^91, (1) 



where the sum is over the 3A quarks, we have for its square, 



The magnitude of any single term is 

?1 = ^/. (3) 

For i ^ j and 3^ quarks, one encounters the same color (repulsive) and different color (attractive) 
interactions 

^ ^ [ 1^2, 3^(^-l)/2 times 

*"^^ = 1 2 2 .2 ^^) 

I -|s''',3A'' times 

respectively. This leads to 



^2_4 2,o, , ,M{A-l) 1. 



g^ = -g\3A + 2(^ ^ - -3A^)) = 0, (5) 

as it should for an overall color singlet. The average repulsive coupling between quarks with the 
same color is 

and the average attractive coupling for different colored quarks in a baryon is 



^_M^V_?„^, = r^ 
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\ 2 / ^ ^ 

The TF quark model replaces the sum over particle number in particle interaction models with 
a sum over the density of state particle properties. The interaction strengths are taken from the 
classical theory, but it is necessary to weight these by the probabilities of the various interactions 
in the color sector, which we assume are flavor independent. We need a connection between 
particle number and probability. The natural assumption is that these interaction probabilities 
are proportional to the number of particle interaction terms. Thus, we assume that the average 
repulsive (same color) or attractive (different color) interaction couplings, Eqs.(|6]) and ([7|, equal 
to ^g'^Yli-^ii ^^"i ~^9'^Yli=/=j^ij^ respectively. In addition, we assume that all colors in a color 
singlet contribute equally to these sums. We therefore have 

^^^ = 3(3^^)' (^^ 



for self-color interactions and 

for different colors. These sum to one as they should: 

5]P,, = 1. (10) 

Using these, we can now construct the nonrelativistic interaction energy, E, from the kinetic and 
color Coulomb parts of the interaction to form the nonrelativistic system energy. This will be done 
in Section III. The overall "residual" coupling, the sum of dm) and ([7]), is 

3A-1 J 3^ 3A-V ^ ' 

which is just the negative of the charge of a single quark divided by the number of remaining 
quarks. 

III. SYSTEM ENERGY AND EQUATIONS 

Assuming flavor and color number, nl{pp), and non-interacting quark kinetic energy, 81{pf), 
densities, we have that these are related to the Fermi momenta, {pf)1 by 

and 

where the "I" superscript stands for flavor and the "i" subscript stands for color. According to 
the assumptions of the TF model, all quantities are derived from the particle densities, n[. Thus, 
according to Eqs.(12) and (13), the kinetic energy density is proportional to (nj)^'^. The model 
assumes that this relationship holds at each point in space, r, and the ground state of the system is 
obtained by minimizing the energy functional. I assume the flavor-summed spatial normalization 
is given by 

^ I cPrnl{r) = A. (14) 

The number of quarks with flavor / is designated N . The color-summed number is therefore 

Y,Jd'rni{r)=N', (15) 



and total number is 

J2^^ = 3^- (16) 

For convenience, I will also introduce the single-particle normalized density 

This form of the density will be helpful in correctly normalizing the TF quark-quark interaction 
energy when continuum sources are used. 

I use the the color interaction probabilities from Section II to construct the nonrelativistic 
system energy, E. I use the notation "i < j" on the color sums below to avoid double counting. The 
minus signs and factors of 1/2 in front of the color sums involving Pij are building in the interaction 
strengths from Eqs.dm) and M. I double count on the flavor sum, I ^ J, and compensate with a 
factor of two. Using the single particle densities, nl{r), I normalize to the number of terms one 
has from the discrete form of the same- flavor interaction, N^{N^ — l)/2, and the number of terms 
in the different flavor interaction, N N . (I will check on the energy normalizations in the next 
section.) One has 

E = T + U =y / dV> *^ " 



i,I 



4 2^^N^(N^ -1) r-"- /-r-™. ^3y,^3^/ / 

+ 3^ E ^—2 / / uFTf^ E ^««"f W"f (- ) 

-^E^^^-^'w^'(^')) 

i<j ) 

+ 3^ E -^T- / / F^Tl \Y,Punl{r)ni{r') 

-^E^*^-^'W^/(^'))- (18) 

I have assumed that the radius of the objects introduced are finite, and indicate this with the 
^^rmax" notation on the spatial integrals. I now do two things: switch to normalization nj and 
assume equal Fermi color momenta, n = n{ = 71-2 = n^, for each /. Doing the color sums, we have 

3(37r2/i3n^(r))5/3 



E = Y. d\ 



X l^^ s-^N^-l r™- T'""" 3 ,3 /n^(r)n^(r') 



2(3>1-1)^^ N^ 



X 3^ 









— T 

2{3A - 1) f-^. 

J^A^(3 / '""^ d^rn\r)-N^], 



(19) 

where I have introduced Lagrange multiphers, A^, associated with the constraint 

I ^^^ d^rn\r)=N^ /?,. (20) 

I now do the variation of the density, 6n^{r), in E. (The steps of setting n[ = n2 = nf^ and doing 
the variation are interchangeable.) The result may be written 



{p'f? _ _^i + 3 X |g^^ I N^ /^^""V/^VI + ^ 
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(3A-1) I N^ 
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(21) 



with pp given in terms of n^ by Eq.(12). Notice the residual coupling, Eq.(ll), as an overall factor 
in front of the potential term on the right hand side. I introduce 



fir) 
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2 X |a, 



-(37r2n^(r))2/3, 



for the TF spatial functions, where 



is the reduced Compton wavelength and 
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OLx 



(22) 



(23) 



(24) 



defines the strong coupling constant. I choose the scale "a" to be associated with the lightest quark 
mass, which I have designated as m}. Since the masses, rn, , are not necessarily equal, I will choose 
to order, for example, m} < m? < m^. Assuming spherical symmetry 



d3/i^A!LL=47r 

\r — r'\ 



dr'r'^ — ^^ + / dr'r'^ ^ ' 



(25) 



one now has the quark model TF integral equations. 
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(26) 



where a^ = m}- jvn} . It is convenient to introduce a diniensionless distance, x, given by 



r = Rx, 



(27) 



re the physical distance R is given by 
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now have 
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The first derivative of this is 

idf^x) X^R A 
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dx 
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3^' ^ 3 



and the second derivative is 
id^f'ix) 



a 



dx"^ 



A 1 
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The normalization integral, Eq.(20), reads 



nax ]\f 

dxVx{f{x)f'^ = —. 

"J^ 



(28) 



(29) 



(30) 



(31) 



(32) 



IV. NONRELATIVISTIC iV/-FLAVOR MODEL 



I will solve the TF inegral equations, (26), for the convenient case of Nf mass-degenerate quark 
flavors with equal numbers, N , even though it may not be realistic. One would expect increased 
stability for larger Nf from Fermi statistics and this will be examined below. We have 

. 3A 



N' 



Nf 



(33) 
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Defining n(r) = n^{r), f{x) = f^{x), and A = A^, where from Eq.(22| 

3 



n{r) 



1 / 2 X |as 



37r2 



we have 



d rn{r) 



SttAJ \xJ 



A 



3/2 



Nf 



In terms of the TF function, /(x), one has 

/ dXy/^{f{x)f/^ 

Jo 

The A'^y-degenerate TF integral equation becomes 

XRx 



Nf 



fix) 



3^ 



+ iV 



/ 



, / ,2 1 fW) 

ax X ' 

x' 



3/2 



+ X 



"% / / [ fix') 

ax X 



The first and second derivatives are 

df{x) XR 



dx 



-^^ + Nf dxx [ — — 

tg^ Jx \ x' 



3/2 



and 



d'fjx) 
dx"^ 



-N 



f- 



(/(x))3/2 



Using Eq.(39|, the normalization integral may be integrated to 

df 



dxV^(/(x))3/2 



which gives the endpoint differential condition 



dx 



1 



Nf \ dx 



XR 

3^ 



where A is the Lagrange multiplier. Notice that the scaling substitutions, 



f^f,x 



N 



2/3'' 
/ 



Nfx 



completely converts Eqs.(36)-(41 ) into the equivalent single flavor system. 



3/2' 



(34) 



(35) 



(36) 



(37) 



(38) 



(39) 



(40) 



(41) 



(42) 



The equations of motion for this model may be solved numerically. The TF function obtained 
for Nf = 1 is shown in Fig. [T] Since there is no central Coulombic source, the value of the function 
is zero at the origin, as can be understood from Eq.lpTJ). It extends out to Xmax = 5.0965. The 
particle density function, which is proportional to {f{x)/x)^'^ (see Eq.(22)) is shown in Fig. [2j also 
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f(x) 




FIG. 1: The dimcnsionlcss TF spatial function, f{x), related to the particle density by Eq.(22l, for the non- 
relativistic Nf = 1 flavor model. The dimensionless distance, x, is related to the physical radial coordinate, 



r, by Eqs.(27l and (28 1 



( f(x) /x 



1.5 




FIG. 2: The (unnormalized) particle density function, {f{x)/x)^^'^ (see Eq.(22)) for the nonrelativistic 
Nf = 1 flavor model. 



for the Nf = 1 case. It has a smooth profile, even at the surface. Remember that the general Nf 



case can be recovered from these results from the inverse of the (42) substitutions. 



I now wish to evaluate the kinetic and potential energies of this model. We will find that these 
can also be related to the endpoint derivative, {df/dx)x,^^^- For the kinetic energy, T, we start 
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with 

T = Nj 
In terms of the TF function, f{x), one has 



3 3{3TT^h^n{r)f/^ 



Wn'^h^ 
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4 y a 



or, using Eq.(39) 
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57r 



24 f^TTAV^^ y^ X |a. 
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Finally, using Eq.(41) we have the simple form 

1/3 4^2 y 4 



^_24 /37r^y/^ |ff^x|«, _9 /|gM 

/Vr y 4 y CLXmax ' \ fmax 



The potential energy, [/, is given by 
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3x|ff^^/ f f^3^^3^ Mr)n{r') 
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(43) 



(44) 



(45) 



(46) 



(47) 



This is probably the appropriate time to explain the philosophy behind the choice of energy/particle 
density normalization in the model. The discrete form of the Coulomb interaction between 3A 
colored objects is 

1 



U 



discrete 






1)5: 



2 X (3A-1) ^^ \ri -rn 



(48) 



We see in Eq.(48| the residual interaction coupling, discussed in Section II, times the Coulomb 

1^ 



interaction of charges with normalization g ^ = ^g^. The number of terms is 3yl(3^ — 1), given the 



double counting in the sum. We may write Eq.(48) as 



^n^ X 3/1 

Tjdiscrete _ 3^ < U > 



(49) 



where < U > represents the average interaction potential energy term. On the other hand, if the 
single particle norm h = n-^ [J d^rh{r) = 1) is used in (47), we have 



U 



1^2x3^ /■ /■ ^3. ^3/^(0^(0 



dV d-'r 






The comparison yields 



<[/>< 



— — —n[r)n[r ). 
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FIG. 3: The dimcnsionless potential, $ = —f{x)/x, from Eq.(53 1, associated with the nonrelativistic Nj = 1 
flavor model. 



This is a natural association, which I consider the appropriate way to to set the overall energy 



normalization in the model. Of course, (48) does not contain self-energy contributions whereas 



(47) does, so the comparison can not be exact. 



I use Eq.(21) to define the potential, V{r), in this model: 



V{r) = A 



A 



3 n{r) 



d-'r 






Evaluation gives 



V{r) 



R X 



(52) 



(53) 



where I have used the definition of the Lagrange multiplier, Eq.(41). Its role is to remove the 



discontinuity in the potential at the surface. The dimensionless form of this potential is displayed 
in Fig. |3j It shows an attractive collective potential with no central Coulombic spike. 
Relating the potential energy to the TF function, /(x), we have 
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Using Eq.(4T| 


, we have 
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as we expect, 


and the energy 


of the system is just 
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The system energy is proportional to A^'^, which means it is always energetically favorable to 
populate n (n = 1,2,3...) A = 1 states rather than a single A = n state. Of course all these 
states would be occupied at nonzero temperature. The size of these loosely bound systems grows 



like ~ A^i^, according to Eqs.(27) and (28). This is in contrast to the usual increase ~ A^'^ for 



nuclear systems and is a result of the "residual" nature of the coupling as discussed in Section II. 



Also note that although the number of degenerate flavors, Nf, does not appear explicitly in (58), 
the scaling substitution, (42), assures that the Nf flavor energy, E^f , is related to the single flavor 



energy, E^^, by E i = N J E^. Since all energies are negative, the A^/-degenerate system is more 
strongly bound than the single flavor system for the same value of A, as it should be. 

Although interesting in themselves, it is clear that the states being formed do not constitute 
normal hadronic matter. As discussed in the Introduction, I interpret these states as loosely bound 
hadronic states in the zero temperature quark-gluon plasma phase. As a model of the confined 
phase, I next consider what happens to these solutions when a confining vacuum energy density 
term is added, as in the MIT quark model. 

V. NONRELATIVISTIC iV/-FLAVOR MODEL WITH VOLUME PRESSURE 

I now introduce the volume energy, 

Evol = Y^'iaxB, (59) 

where B is the MIT bag constant. The total energy is of course 

Etot = E^oi + T + U. (60) 

We will again solve for the case of Nf mass-degenerate quark flavors with equal numbers, 
A^^. One may anticipate that a discontinuity will be produced by the introduction of an external 
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FIG. 4: The nonrelativistic TF functions, f{x), as a function of x for Nf — 1. Reading from bottom to top, 
this shows the /{x^ax) = 0; ^: i ^^^d 1 functions, corresponding to x„iax = 5.097,2.741, 1.945 and 1.691. 



pressure, so that f{xmax) 7^ 0. Given this, the normahzation integral is now (use the differential 



equation (39) in the integral and do an integration by parts) 

dXy/x{f{x)f''^ 



Nf V dx 



(61) 



where of course the integral is required to have the same value as in Eq.(36). This leads to the 



modified boundary condition (also use Eq.(38) to relate this to A at the boundary) 

7_r 



dx 



XR 

3^ 



(62) 



These solutions can be characterized by the value of the TF function at the endpoint, f{xmax)- 
These are displayed in Fig. El for the Nj = 1 case. Note how the endpoint value of the TF function 
rises with decreased x^ax- This is shown in Fig. |5] 

Considering that there is a pressure-related discontinuity at the surface, the kinetic energy in 
this model is now given by 



24 f^y^^g^x las 



p^-rm -^^/- 
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UsingEq.(62| 
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Likewise, for the potential energy, we have 
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FIG. 5: f{xmax) as a function of Xmax for the Nf = 1 case. 



Again using ( 62 ) 
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and therefore 



E = T + U 
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4 /37ryl\ 



vr 



i^^ + ^VS/V^ 
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I wish to minimize the energy with respect to rmax, which, since A is held constant, is the 
same as minimizing with respect to the dimensionless radius, Xmax- In taking the derivative of the 
total energy, one has to be careful about the meaning of derivatives. One is taking the derivatives 
of f{Xr 



with the normalization integral, Eq.(36), being held constant. These special types of 



derivatives will be called "endpoint" derivatives to distinguish them from the usual derivatives 
applied internally to the TF functions. They will be designated with a subscript, "e". Using this 
notation, the pressure matching condition is 

dEtot 



dXr 



0. 



(68) 



The numerical relationship between f{xmax) and Xmax is given graphically in Fig.l5](A''j = 1 case). 
Several limits of this figure may be easily understood. At the final point, Xmax = Xf {xj = 5.0965 
for f{xf) = 0), one may show that {df{xmax)/dxmax)e\xrr^aa,=xf = -l/a^/- The value of the endpoint 
derivative of the function in Fig. [5] is displayed in Fig. [6] for 2 < x^ax < Xf. Note from this figure 
that there is a small maximum near the final point, x/, leading to a small initial decrease in the 
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-0.2 



-0.25 



-0.35 




FIG. 6: The endpoint derivative of the function from Fig.lslin the region 2 < Xmax < Xf. Note the maxinium 
in the function near Xmn-r — 4.5. 



slope. It is otherwise monotonically decreasing. In the other hmit, x„ 
is obtained as A — > oo, one may show that 



fixr, 



U 



2/3 



1 



0, which we will see 



(69) 



Taking the derivatives of the various terms and setting them equal to zero to balance the 
pressure gives 



^2 (Ba^/hc\ (?,ttAV'''^ 



ai 




7x \ dx 



35 ^/x V \dx J ^ 



(70) 



where again we encounter endpoint derivatives. I will call this the "matching equation" and the 
right-hand side the "matching function" . It is displayed in Fig. 17] for 4 < Xmax < Xf . Note that 
the function is initially negative in the vicinity of x^ax = Xf. This means there are no stable 
solutions in this region. This instability is a consequence of the small initial decrease in slope seen 
in Fig. |6| The negative excluded region decreases as Nf increases. It excludes the approximate 
regions 4.83 < Xmax < xj for Nj = 1, 5.03 < Xmax < Xf for Nj = 2 and 5.07 < Xmax < Xf 
for Nf = 3. Although it continues to shrink, the excluded region seems to persist for any value 
of Nf. Truly heavy quark systems, such as for charmed or bottom, would have small "a" values 
that would force them to have Xmax values near the edge of this region. Finally, notice that the 



left-hand side of Eq.(70) is obviously a very sensitive function of both a and a^. 
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FIG. 7: The "matching function" of Eq.(70l a function of Xmax for Nf = 1 in the region 4 < Xmax < Xf. 



Using Eq.(70) in (59) allows one to formally eliminate the bag constant B, and gives the total 
energy as 



E, 



tot 



SttAV^^ (Ig^ X las^ 



a 



8/-1 



2 /# 
7 Ux 



)y-^v-.f'if 



X 



df 



(71) 



7 X 7 \ dx J „ 21 V \dx ^ 

Note that Etot = is not the reference energy for a system being bound, since there are no free 
particle states. One can only compare the energies of one set of bound states with another in the 
confined phase. 



The energies as a function of A from Eq.(71) are complicated and they have not yet been 
explored. However, in the limit of A — > cxd, one can show from Eqs.(69) and (70) that the system 
shrinks, Xmax ~ A~^'^, leading to Etot ~ A. However, this limit is unphysical because the system 
is becoming relativistic as it is being confined to a smaller volume. This leads us to our next set 
of considerations. 



VI. RELATIVISTIC MODELS 



Can the present model be made relativistic? Although the expression for the relativistic energy 



density, {£r)1, is more complicated, the equation replacing Eq.(21) above is clear: 

'iV^-1 



{p^pc)'^ + {mic^Y ~ ""^-fC 



J 3 X I52 



^R 



(3A-1) 



Ni 



,3 n\r') 



d-'r 
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n.ax /y^J ( fy^\ 






(72) 



with pp again given in terms of n by Eq.( 12 ). This reduces to the previous nonrelativistic equation 
in the low density hmit. Instead, I will concentrate on showing the existence and properties of the 
extreme relativistic or massless version of this equation. This is given by 



T .T 3 X 1^2 / ^/ _ 1 



,3 n\r' 



(fr 















(73) 



I will continue to assume A'^^ degenerate flavors with equal numbers, A^ , in this section. Quantities 
such as p^p, A|j, and n^{r) now lose their flavor superscript. I will introduce the notation w{r) for 
the degenerate massless TF function. It is dimensionless and is related to the previous definition 
of /(r), assuming m 7^ 0, by 

1/2 



w{r) = ^(3vr2n(r))^/^ 



2mc f{r)r 



la. 



3^ 



(74) 



Eq. ( 35 ) now reads 



rmax ^^ ^ ^ ^^^ I 



where 



, (-(0)^ - ^, 



/3 = 



^{hsYNf 



?,-kA 



(75) 



(76) 



We follow a path similar to the one for the nonrelativistic model. The TF integral equation in 
this case reads 



w{r) 
The first derivative is 



Art 

4^2 



+ /3 



dr' 



ax^^l 



-(u'(r')) +r / -Piiwi'r')) 



— --43+/3/ ;72 W)) 



dr 



^r 



while the second is 



d^w 
dr"^ 



-4. 



Eq.(75) nows gives 



fdw^ 



w — 1 



4^2' 
3^ 



(77) 



(78) 



(79) 



(80) 
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FIG. 8: The relativistic TF wavefunctions as a function of x. Reading from bottom to top, this shows the 



w{rmax) — 0' li i ^^d 1 functions with r„ 



1. 



similar to Eq.(41 ). A version of Eq.(|79[) appears in Ref.[Tl] in the context of a study of the quark- 
gluon plasma with SU{2) or SU{2>) color. In this paper the authors subject the partition function 
to a TF approximation and confine the system of massless quarks with an external, spherical box 
in order to study its thermodynamics. Instead of the single function w{r)^ this development has 
two potentials, called /(r) and a{r)^ for SU{2). If the gauge particles had abelian f/(l) symmetry, 
as is assumed here, the o(r) function would not arise. The /(r) equation would then go over 



to a scaled version of Eq.(79) in the zero temperature, zero density (chemical potential) limit. 
However, the boundary conditions used for the field functions are quite different from those used 
here. Presumably, the equations presented in Ref.[l4, will not contain finite quark matter solutions 
since (3 is replaced with a negative constant, producing repulsion, and /(r) is seen to diverge for 
large distances. Perhaps requiring the plasma to be an overall color singlet, as is done here, would 
change the sign of this coefficient. 

Getting back to Eq.(79), notice that it is invariant under the substitution r — > Sr, where S 



is a scale factor. Thus Eq.(77) formally has a solution, but no scale! Unlike the nonrelativistic 
model which has the Compton length as an intrinsic scale, using w{rmax) = as a boundary 
condition is meaningless here. This will be overcome by introducing the bag constant, which, like 
the nonrelativistic case, will result in a discontinuity in the TF function at the boundary. Formal 
solutions for the TF function w{r) are shown in Fig. ^ where I have set rmax = 1 in arbitrary 
units. 
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The ultra-relativistic kinetic energy is given by 
or, when related to w{r), we have 

VT 3 ./o 



""" 3 3c(37r2;i3n(r))4/3 



4^2^13 



'\aX ^rp 



{w{r))\ 



(81) 



(82) 



Evaluation, using Eq.(79), relates this integral to endpoint values and derivatives of the w function: 

9.4 






(83) 



Using Eq. ( 80 ) then gives 



--yh'^- 



1-W + Pw^/2 



(84) 



It is interesting to examine the relativistic potential, VR(r). Similar to Eq.(52) above, one has 



VR{r) = \r 



A 



d-'r' 






This results in 



VR{r) 



4 ^wir) 
3 r 



(85) 



(86) 



The relativistic Lagrange multiplier, Ajj, plays the same role here as in the nonrelativistic model. 
The density function and the potential are displayed for the w{rraax) = version of the model in 
Figs. |9]and 10 setting r^aax = 1 in arbitrary units. One again sees smooth density and potential 
profiles with no central Coulombic charge. Of course, the 'w{rmax) / models inevitably involve a 
discontinuity in these quantities at the surface. 

The potential energy is (formally the same as Eq.(|47|)) 
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3 X y^Nj 
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rfVcivl"^"!"!^ 






(87) 



When related to the TF w function, this becomes 
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Using techniques similar to the nonrelativistic evaluations, one obtains 
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(88) 



(89) 
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FIG. 9: The (unnormalized) relativistic TF density, proportional to {w{r)/r)^ from Eq.(74l, for the scaleless 
w{rjnax) = model with r„iax = 1- 




FIG. 10: The relativistic TF potential, ^r = —w{r)/r, from Eq.(86l, for the scaleless 'w{rmax) — model 



with r„^ax = 1- 



Finally, using Eq. ( 80 ) , this is equivalent to 

9,4 



U 



,qi)A 



1-W + pw^/3 



(90) 



The total energy of the unconfined system is given by 

chNf A 



E = T + U 



vr '3 



a. 



w 



(91) 
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FIG. 11: /3 as a function of w{rjnax)- 



The situation here is intrinsically different from the nonrelativistic model. The solution to 



Eq.(79) with endpoint boundary conditions given by Eq.(80), for arbitrary r^ax^ shows that there 



is a one-to-one relationship between (3 and w{rmax)- This connection occurs of course because one 



is keeping the normalization integral, Eq.(75), constant as the surface discontinuity is varied. One 



finds that (if = 4.073 corresponds to w{rmax) = 0. This is illustrated in Fig. 11 where (i is plotted 
as a function of w{rmax)- Notice that the function in this figure has zero slope at w{rmax) = 0. In 
the other limit, Wmax — ^ 0, which is obtained as A — > cx), one can show that 



/3 



w 



(92) 



For the confined system we define the total energy as in Eqs.(59) and (60). I again seek the 



minimum energy solution by varying the radius, rmax while keeping A, and therefore w{rraax), 
constant. This gives 



dE, 



tot 



<9r„ 



47ri?r^ 
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(93) 



max / e 

By requiring this to be stationary one obtains a linear relation between rmax and w(r„ 



-.a. 



Nf ch 



1/4 



'^ymax)' 



(94) 



3'" V47r2 B 

Notice that the A value does not appear explicitly here, but is implicitly involved because of the 
connection between w{rmax) and p. The radius formula is no longer a matching condition to be 
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solved, as in the nonrelativistic case, but simply a pressure matching result for the radius. After 
eliminating the bag constant, the total energy can be written as 

Etot = ^{-^a.) ^-j . (95) 

^ '^ ' max 

Energies here are intrinsically positive. The relationship between (3 and w{rmax) implies a nontrivial 



A dependence in Eq.(95) from the relationship between w{rmax) and (3. Again, I have not yet 



explored the A dependence of the energy in general. However, in the the A — > oo limit, Eq.(92) 
gives Efot ~ Z?""*^) which shows that the total energy is proportional to A and inversely proportional 

to Nf. 

VII. COLOR-FLAVOR LOCKING CONSIDERATIONS 

One of the interesting possible applications of the present model would be to study so-called 
quark strangeletsjll [2], hypothetical finite lumps of hadronic matter made stable by an optimal 
mixture of up, down and strange quarks. Much work on such systems has already been done, 
especially in bag mo dels [T7]. The question of stability is critical, and comparison of the TF quark 
model with previous results should provide a consistency check. 

In a very interesting set of papers, the authors in Ref.[r8| (see also the review articles [T9]) have 
given evidence of a scenario in which the ground state of cold, low density quark matter is described 
by a mechanism called color-flavor locking (CFL). It is conceivable that for a large number of quarks 
this could produce finite quark matter[20]. CFL is due to an attractive instanton or single gluon 
exchange interaction between two quarks with different flavors and colors, and results in Cooper 
pairing and a diquark condensate. It also leads to an equal density of up, down and strange quarks, 
at least in the zero quark mass limit. Moreover, the spontaneous symmetry breaking involved leads 
to massive gluons via the Higgs mechanism. I will examine some aspects of the CFL modifications 
on the massless quark solutions of the last section, although these changes could also be made in 
the context of the nonrelativistic model of Sections IV and V. 

Let us define /i = mgc/h, where rug is the gluon mass. I assume A'^^ = 3 degenerate flavors 
and introduce a Cooper pairing gap. A, in the momentum spectrum of the massless quarks. The 
formation of Cooper pairs leads to an approximate energy density [21], 



Opair - - ^^3^2 ' ^^") 

where A is the energy gap. In the TF model, such a term becomes a function of the color-flavor 



density function, n{r), using Eq.(12), and can be incorporated in the energy functional. 
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FIG. 12: TF massive gluon functions, w{r), for (a)A* = 10 (b)A* = 1 (c)A* = 0.1 (d)A* = 0.01. A* in units 
of hc^ and distances measured in units of /i^^. 



The modifications to Eq.(77) of the last section are now straightforward. We find 
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(97) 



Of course, the r dependent quantities may be taken outside the integrals in (97). One does not 



expect a singularity on the left-hand side from the pairing term due to the linear vanishing of w{r) 
at the origin, and solutions with w{rmax) = are avoided when the bag constant is employed. As 
in the last section, I will bypass a general consideration of this equation in favor of a special case. 
We will concentrate on looking at the massive gluon modifications and set A = in the following. 



The first derivative of Eq.(97) for A = gives 
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At rm.ax one has 
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The second derivative is then 
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FIG. 13: /3 as a function of A*, measured in units of hc^, for the massive gluon model. The dashed horizonal 
hne gives the asymptote /3f — 4.0733. 
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FIG. 14: A log-log plot of r^ax, in units oi n ^, as a function of A*, measured in units of hciJ, for the 
massive gluon model. 



Notice the extra terms involving //, which now gives this massless quark equation a scale. 

I only consider the the zero pressure case which has w^rmax) = 0. The solutions of these 



equations using Eq.(99) as a boundary condition for various values of A* are illustrated in Fig. 12 



Note that in the limit of large A* the TF functions shrink in size and begin to look suspiciously 
like the TF w{rmax) = function in Fig. pj This is not a coincidence. In fact, one can show that 



Eq.(lOO) goes over to the simpler Eq.(^79p in this limit. On the other hand, for small A*, the TF 
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function grows in size and begins to develop a longer tail. The relationship between A* and /5 is 



given in a log-linear plot in Fig. 13 Note that as A* — > cxd, we recover the /?/ = 4.073 limit 



for the unconfined massless case considered in the previous section. The relationship between A* 



and rjTiax is shown in a log-log plot in Fig. 14 From the solutions to Eq.(lOO) one can show that 
fmax — > 1/A* in the this limit. This is evident from the log-log plot of Fig. 14 which gives a slope 
of -1 in the A* — > cx) limit. 

The energy functions can no longer be formally integrated due to the gluon mass term and I 
will not attempt to study the A dependence of the energy here. 

VIII. SUMMARY AND COMMENTS 

I have introduced a semi-classical TF model of quark interactions in finite quark matter. A MIT 
like vacuum energy density is used to model confinement, but the (massive) model itself provides 
effective confinement |22] from the attractive Coulomb interactions, at least in the absence of CFL 
locking. I have presented nonrelativistic and relativistic forms of the model and shown that it can 
also be adapted to the CFL scenario. 

I have computed, compiled, and derived a large number of results and relationships, including 
the TF functions, densities, potentials and other quantities which characterize the various models 
introduced. However, it must be admitted, in the presentation of topics for investigation I have 
picked the "easy" problems to solve. I have looked at the nonrelativistic case without examining the 
more difficult coupled equations for unequal quark masses. I have also solved the ultra-relativistic 
limit after bypassing the much harder massive relativistic equation. Finally, I looked at the CFL 
modifications for massless quarks without attempting to solve the gap equation for Cooper pairing 
or incorporating the bag constant. Clearly, one can not make contact with phenomenology without 
a significant amount of further work on the harder problems. However, the point of this work is 
simply to show the range of problems which can be addressed and the relative ease of applicabihty 
of the TF model for quarks. Future effort will hopefully connect with the important questions of 
relative stability and particle phenomenology. Although the TF model assumes a large number 
of particles, the amazingly accurate binding energies found in the atomic case for low atomic 
number [12] gives us hope that the results obtained from this model will also be reliable for small 
A. 

A word about parameters is in order here. I do not expect to fit hadron masses and other 
measured quantities in order to set the values of parameters in the TF model, which are a^. 
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B, and the various quark masses. That is the domain of detailed spectral models. Instead, 
phenomenologically relevant values of these parameters will be used and explored. Hopefully, many 
of the relations, for example, hadron mass or binding energy ratios, will be relatively insensitive to 
these parameters. We will have to see. 

Among the many physical applications I have neglected two deserve special mention. One of 
these is the heavy quark case. In the limit of heavy masses quarks will become nondynamical. 
This can be modeled by the introduction of central Coulombic color sources, very similar to atomic 
nuclei, in the midst of dynamical light quark distributions. Their effect on the stability character- 
istics of the system could be dramatic and such systems are relatively unexplored. In addition, the 
considerations of Section II can easily be generalized to systems with antiquarks. This will involve 
introducing new color couplings in the Coulomb interaction between quarks and antiquarks as well 
as antiquark antiquark interactions in the system energy, E. Of course, any possible annihilation 
diagrams will have to be neglected. Both of these topics will allow a new class of objects to be 
investigated in this model. Now that the foundations of the model are established, further work 
can concentrate on fleshing out some of these neglected topics. 
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